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Abstract 

This talk describes the physics input of QCDINS, a Monte Carlo 
event generator for QCD-instanton induced scattering processes in 
deep-inelastic scattering. 
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Hard scattering processes in strong interactions are snccessfnlly described 
by pertnrbative QCD. However, pertnrbation theory does not exhaust all pos¬ 
sible hard scattering processes: Instantons , non-perturbative fluctuations 
of the gluon fields, induce hard processes which are absent in perturbative 
QCD. Deep-inelastic scattering (DIS) at HERA offers a unique window to 
detect these processes through their characteristic multi-particle final-state 
signature P|. A Monte-Carlo generator for instanton-induced events in DIS, 
QCDINS, interfaced with HERWIG, has been developed 0 which enabled 
the HI Collaboration to place first experimental upper limits on the cross- 
section [^] and which allows for the elaboration of dedicated search strategies 
1^. The purpose of this talk is to outline the basic physics input of QCDINS 
and its built-in features, characteristic for the underlying instanton mech¬ 
anism. Further details about the theoretical background p and ongoing 
experimental searches appear elsewhere in these proceedings. 

In order to appreciate the notion of instanton-induced scattering pro¬ 
cesses, let us recall that scattering amplitudes are derived via analytic con¬ 
tinuation and LSZ reduction from Euclidean Green’s functions, which in turn 
can be represented by a path integral, 

^ J[dA][d^/J][d1/J]Af,{xl) ...ipixi) ...ij{xn)exp{-S[A,ij,i/j]}. ( 1 ) 

The perturbative scattering amplitudes are obtained from an expansion of 
Eq. (|l|) about the perturbative-vacuum solution, i.e. vanishing gluon fields, 
A^^l = 0, and vanishing quark fields, = 0, with vanishing Eu¬ 

clidean action = 0. This expansion can be summarized by the familiar 
Feynman rules which construct the perturbative scattering amplitudes in 
terms of propagators and vertices as a power-series in the strong coupling 
(see Fig. (left)). 

The instanton (x) is a non-trivial solution of the Euclidean gluon-field 
equations and thus a non-trivial local minimum of the Euclidean action with 
AH) = 271/as- Instanton-induced scattering amplitudes, being derived from 
an expansion of Eq. (|^ about the instanton, can be constructed according 
to modified Feynman rules, which involve, in addition to the propagators 
(in the /-background) and vertices (in the /-background) also the classical 
fields (see Fig. |l| (right)). Instanton-induced scattering amplitudes are always 
exponentially suppressed at weak coupling, oc exp{—27r/Q;s}. 

In QGD with massless quarks, usual perturbation theory and instan¬ 
ton perturbation theory describe two distinct classes of processes: In usual 
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perturbation theory, quarks are coupled to gauge fields only via vector cou¬ 
plings, which conserve chirality (Q 5 ). Thus quarks always appear in pairs 
with net chirality zero (see Fig. (left)) and chirality is conserved to all or¬ 
ders. In instanton-perturbation theory, on the other hand, only scattering 
amplitudes for processes which violate chirality by AQ 5 = 2nf units receive 
non-vanishing contributions p. This is due to the fact that the classical 
quark solutions always appear in pairs with net chirality two (see Fig. ^ 
(right)). 



Figure 1: Tree amplitudes for 'y*g scattering processes, in which all light 
quark flavors are produced. Left: Usual perturbation theory. Right: 
Instanton-perturbation theory. Lines ending with a blob stand for classi¬ 
cal right-handed quark (t/’^^) and anti-quark solutions; line with a 

blob in the middle denotes the quark propagator in the /-background; curly 
lines ending with a blob stand for classical instanton gluon fields 


The Monte-Carlo simulation of /-induced events proceeds in three steps. 
First, quasi-free partons are produced by QCDINS with the distributions 
prescribed by the hard process matrix elements. Next, these primary partons 
give rise to parton showers, as described by HERWIG. Finally, the showers 
are converted into hadrons, again within HERWIG. 

The momentum-space structure of /-induced hard processes in ^*9 scat¬ 
tering is shown in Fig. |] (left). As is already suggested by the form of 
the leading-order matrix elements in Fig. |l| (right), the amplitudes factor¬ 
ize into a product of an effective 'y*qq* vertex 0 , denoted by a blob in 
Fig. I (left), times matrix elements for /-induced partonic subprocesses, 
q*g —> {2nf — 1) q + Ugg, denoted by a blob with an index “/”. It turns out 
that the most important kinematical variables determining the final state 
properties of /-induced events are the virtuality of the off-shell quark q*, 
Q'^ = —> 0, and the Bjorken scaling variable of the q*g subprocess, 
x' = Q'"^/{2p ■ q'). Therefore let us discuss their distributions first. 
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Figure 2: Left: Momentum-space structure of /-induced partonic processes in 
7 * 5 ^ scattering. Middle: The holy-grail function/^(a;'). Right: x'distribution 
from QCDINS. 


To this end, we square the amplitudes in Fig. ^ (left) and sum over 
all hnal state partons to construct the /-induced hard-scattering partonic 
total cross-section or structure function. It can be shown [Q that the initial 
state collinear singularities, occuring when Q'^ —>■ 0, can be consistently 
absorbed into the parton distribution functions fk, such that the /-induced 
contribution to the nucleon-structure function is obtained in the familiar 
convolution form, 

k .Td; 


dx 
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fk 
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where denotes the renormalization (factorization) scale. The I- 

contribution to the dominating gluon-coefficient function in turn, has 
for large the anticipated momentum-space structure Q, 



The essential instanton dynamics and in particular most of the dependence 


4 












on x' and Q' is encoded in the J-subprocess total cross section, 
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where the functions S(a;') and F{x') (see Fig. ^ (middle)) are known Bl 
for x' > x':„ ~ 0.2 — 0.3. 
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We see from Eq. (^ that the summation over the gluon emission has mod- 
ihed the exponential suppression factor by the so-called “holy-grail” function 
F{x'). This exponentiation is due to the fact that every gluon emission brings 
in a factor of ~ 1/^/al in the exclusive amplitudes (see Fig. E (right)), 


such that 0 


a, 


(I) 

q’-g-^rig 


cx (l/ng!)(l/as)”9 exp{-47r/a4- 


The holy-grail function F{x') is continuously decreasing from 1 at a:' = 1 
(low q*g c.m. energy) to 1/2 at x' ~ 0.2 (see Fig. (middle)). The total cross- 
section is correpondingly exponentially growing with decreasing x'. This is 
clearly seen in the x' distribution taken from the Monte-Carlo simulation (see 
Fig. I (right)). 

Thanks to the chosen renormalization scale, g.{Q') = Q'as(/r(Q'))/(47r), 
the J-subprocess cross-section (||), as a function of the q* virtuality Q', has 
a peak structure which is clearly reflected by the Q' distribution in Fig. ^ 
(left). The peak is at around 5 GeV, which is gratifying since it means that 
even without a lower Q' cut[| most of the events generated are hard enough 
in order to justify instanton-perturbation theory 0. 

In Fig. ^ (right) we present the resulting /-induced total eP cross-section 
for HERA, subject to the following cuts: i) xbj > min, VB] > 0.1 in 
Eq. (^); ii) x' > xF^ in the integration of Eq. (|^). No cut on Q' has been 
imposed. It is important to note that lower limits on x' and Q' can be 
enforced experimentally by cuts on hnal state momenta0. The points in 
Fig. I (right) have been taken from the Monte-Carlo simulation (QCDINS 
1.4.1). As a check of the Monte-Carlo, we have also analytically integrated 
Eqs. (^ and (|). The resulting cross-section (lines in Fig. ^ (right)) nicely 
agrees with the Monte-Carlo result. We refrain from going to even smaller 
XBj, say 10“^, since in this case = xbj Vb] S would be only of order 1 GeV^. 
It is not clear whether in this case the corrections to Eq. (|^), which are not 
built into QCDINS, can be neglected. 


minimum Q' can be enforced by requiring the transverse momentum kr of the 
current-quark jet (see Fig. || (left)) to be large. A lower a;'-cut can be implemented, for 
example, by requiring x = + s).^ where s is the 7*5 c.m. energy. 
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Figure 3: Left: Q' distribution. Right: /-induced total eP cross-section 
for HERA (preliminary) with various cuts as indicated. Lines: Analytical 
calculation. Points: Monte-Carlo simulation with QCDINS 1.4.1. 


As a benchmark for searches for /-induced events at HERA H.0. let 
us compare the /-induced HERA cross-section from Fig. ^ (right) with the 
normal DIS HERA cross-section, > 10“^,|/Bj > 0.1) ~ 15 nb. We 

see that the fraction of /-induced events at HERA, 
ranges between 

0.0002 % </(^)< 0.1%, (5) 

in the kinematical range x' > 0.3..0.2, xbj > 10“^, i/bj > 0.1. Note that the 
published upper limits on the fraction of /-induced events placed by the HI 
Collaboration [Q are in the several percent range. An improvement of these 
limits by an order of magnitude is reported in another talk in this working 
group 0. 

Let us turn now to the hnal states of /-induced events in DIS. The current 
quark in Fig. § (left) will give rise, after hadronization, to a current-quark 
jet. The partons from the /-subprocesses, q*g {2nf — 1) q + Ugg (see Fig. g 
(left)), on the other hand, are emitted spherically symmetric in the q*g c.m. 
system (“/-c.m. system”). The gluon multiplicities are generated according 
to a Poisson distribution with mean multiplicity 


{ng{x', ~ 


27r 


as(h(Q')) 


, ,^dF{x') 


( 6 ) 


which is of the order of two for x' ~ 0.2, Q' ~ 5 GeV (see Fig. ^ (left)). The 
total mean parton multiplicity is large, of the order of ten. After hadroniza- 
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Figure 4: Left: Gluon multiplicity distribution of /-induced events. Right; 
Lego plot of a typical /-induced event in the HERA-lab system at xbj = 10“^. 


tion we therefore expect from the /-subprocess a hnal state structure remi¬ 
niscent of a decaying hreball: Production of the order of 20 hadrons, always 
containing strange and possibly charmed mesons, concentrated in a “band” 
at hxed pseudorapidity rj in the (r^, azimuth angle 0)-plane. Due to the boost 
from the /-c.m. system to the HERA-lab system, the center of the band is 
shifted in r] away from zero, and its half-width is of order Ar] = 0.9, as typical 
for a spherically symmetric event. The total invariant mass of the /-system, 
\^P = Q'— 1, is expected to be in the 10 GeV range, for x' ~ 0.2, 
Q' ~ 5 GeV. The lego plot of a typical /-induced event shown in Fig. | 
(right) shows that these expectations (current-quark jet; hadronic “band”) 
are actually bourne out from our Monte-Garlo simulation. 
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